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THIN OBLIQUE AIRFOILS AT SUPERSONIC SPEED 

By Robert T. Jones 


SUMMARY 

The well-known methods of thin-airfoil theory have been 
extended to oblique or sweptback airfoils of finite aspect ratio 
moving at supersonic speeds. The cases considered thus far are 
symmetrical airfoils at zero lift having plan forms bounded by 
straight lines. Because of the conical form of the elementary 
flow fields, the results are comparable in simplicity to the 
results of the two-dimensional thin-airfoil theory for subsonic 
speeds. 

In the case of untapered airfoils swept back behind the 
Mach cone the pressure distribution at the center section is 
similar to that given by the Ackeret theory for a straight airfoil. 
With increasing distance from the center section the distribution 
approaches the form given by the subsonic-flow theory. The 
pressure drag is concentrated chiefly at the center section and 
for long wings a slight negative drag may appear on outboard 
sections. 

INTRODUCTION 

In reference 1 it was pointed out that the wave drag of 
an infinite cylindrical airfoil disappears when the airfoil is 
yawed to an angle greater than the Mach angle. This 
observation led to the conclusion that the drag of a finite 
airfoil could be greatly reduced by the use of sufficient 
sweepback. "With such a sweptback wing the wave drag 
would be associated with departures from the ideal two- 
dimensional flow at the root or tip sections and would thus 
be a function of the aspect ratio. The present report extends 
the theory of reference 1 to take account of these effects. 

The treatment is based on the theory of small disturbances 
in a frictionless compressible fluid. The idealized fluid and 
its equations of motion are identical with those employed in 
acoustics in the theory of sound waves of small amplitude. 
The application of the theory is thus limited to bodies having 
thin cross sections so that the velocity of motion imparted 
to the fluid is small relative to the velocity of sound and so 
that the pressure disturbances produced are small relative 
to the ambient pressure. 

The adaptation of the sound-wave theory to the aero- 
dynamics of moving bodies was suggested many years ago 
by Prandtl. The theory was applied by Ackeret (reference 2) 
to thin airfoils moving at supersonic speed. Ackeret’s 
treatment is limited, however, to infinitely long cylindrical 
airfoils moving transversely. The present theory may be 
considered an extension of Ackeret’s theory to take into 

842951 — 80 19 


account wings of finite span and wings having tapered or 
sweptback plan forms. In the case of sweptback plan forms 
the results are markedly different from those obtained by 
the Ackeret theoiy and approach the values indicated in 
references 1 and 3. 

In reference 4 Busemann describes a method for cal- 
culating the supersonic flow over bodies which produce a 
conical pressure field. Busemann shows that the flow 
around cones of circular cross sections as well as the flow 
around the tip of a rectangular lifting surface satisfies 
this condition. The fact that a great variety of three- 
dimensional flows can be constructed by the superposition of 
conical and cylindrical flow fields leads to an essential sim- 
plification of the airfoil theory at supersonic speeds. 

The present treatment differs from Busemann’s in that it is 
further limited to flat bodies, that is, bodies which are thin 
in both longitudinal and transverse sections. This addi- 
tional restriction leads to a much simpler mathematical 
treatment and one which is applicable to a wide variety of 
airfoil shapes. Symmetrical nonlifting bodies are also 
treated' in reference 5 where use is made of integral expres- 
sions corresponding to the velocity potential of plane-source 
distribution. 

SYMBOLS 

V flight velocity 

M Mach number 

x, y, z coordinates 

% point on A r -axis 

Xi limit of integration 

<t> disturbance-velocity potential 

u, v, w disturbance-velocity components 

u t value of u at x t 

u value of u for conjugate arrangement 

p local pressure 

q dynamic pressure 

p density of air 

P„ Q a Legendre functions 
I source-strength factor 

D differential operator 

G D drag coefficient 

t thickness of wing 

m slope of line source (absolute value) 

c chord of wing 
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THE OBLIQUE LINE SOURCE 

The assumptions of small disturbances and a constant 
velocity of sound throughout the fluid lead to the well- 
known linearized equation for the velocity potential 6 (see 
reference 6) 

( 1 — M i ) 4 > xx -\- tf > vv -\-^> tt =0 ( 1 ) 


the other components may be obtained by integrating the 
given cqmponent“to obtain the velocity potential and then 
differentiating the results along the desired directions to ob- 
tain the. desired components. This procedure is especially 
useful in the thin-airfoil theory, where the complete velocity 
field may not be required. 

Adopting the foregoing procedure, one may write 


The analysis is simplified by introducing the coordinates 

*i=* \ 


U o=v 


yi=-jte 2 —l y\ 

Zi= V-Af 2 — 1 si 


( 2 ) 


Dropping the subscripts from the transformed coordinates 
gives 

4>xx—<l>n — 4>it — 0 ( 3 ) 


According to the thin-airfoil theory the pressures on the 
transformed airfoil are given by 


Since u is proportional to the pressure, such a solution cor- 
responds to a point source in the pressure field. The 
solution for an oblique line source may be obtained by 
integrating for the effect of a row of point sources along tho 
line y = mx. It will be shown that such a line source satisfies 
the boundary condition for a thin wedge-shape body. This 
solution, as well as other expressions relating to oblique air- 
foils, can be most conveniently expressed by referring to the 
oblique coordinates 

x'=x—my 


&P_ 0 u 
g_~ 2 V 

2 


(z-» 0) (4) 


(irZ 

and the slope of the airfoil surface ^ is equal to the slope 
of the streamlines near the chord plane; that is, 
dz w 

&TV 


-Icty 

~Y?>Z 


(2-»0) (5) 


The use of the coordinate transformation, equation (2), will 
be understood in the following development. The results 
are therefore applicable directly to a Mach number of V2. 
For an equivalent airfoil at another Mach number the y- 
and 2 -coordinates of the surface will be multiplied by 
\JM 2 — 1 while the pressure coefficients at corresponding 
points will be divided by the quantity M*—l. 

The elementary solution of equation (3) for a point source 
is 

1 . 

This solution is directly related to the subsonic potential 

1 


In the subsonic case the equipotential surfaces are, however, 
ellipsoids, whereas in the supersonic case the equipotential 
surfaces are hyperboloids limited by the Mach cone. (See 
reference 6 for the derivation of these elementary solutions.) 

Because of the linearity of equation (1) a solution may be 
used to denote one of the velocity components rather than 
the velocity potential. The specification of one component 
in this manner actually describes the whole flow field since 


y'—y—mx 

(See fig. 1.) It may be shown that if any function /(r, y, z ) 
is a solution of 

fn fm y«=o 

then /(a/, y', z') is also a solution. In particular, tho point- 
source solution becomes 


1 1 1 

Va •, n —y n —z' t y 1 — m s Vr 1 — y*— 2* 



Figure l.— Oblique coordinates. 
x’-x—my 
y'-y—mx 
2 1 - Vi -in* r 
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Hence the integration for the effect of an inclined line of 
sources may be performed directly along the oblique aj'-asds; 
thus, for m<1.0 


T f { 'i d? 

VO 

x' 

=7 cosh 1 z n (6) 

where is the position of the last source whose Mach cone 
includes the point {x f , if, z') and is given by 

r I =x / -^/y 75 +i 7i 

When m approaches 1.0 the source line approaches coin- 
cidence with the Mach cone, corresponding to a transverse 
velocity component equal to the velocity of sound. 

For values of m greater than 1.0 the integration yields 

u=-Ii . (7) 

Vy /2 +z /2 

It will be seen that in this case 7 is imaginary. 

The vertical velocity near 2=0, which determines the 
shape of the boundary, may be determined by integrating u 
with respect to x and then differentiating the resulting 
velocity potential with respect to z; thus (see appendix), 




where \y' | denotes the absolute magnitude of y'—y—mx. 
The pressure is thus constant along the radial lines 



^ 7 =Constant 


( 10 ) 


= ±x|VT^ 2 (8) 

if z -> 0 and y'< CO. If y'> 0, w—0. There is thus a discon- 
tinuity in the vertical velocity of the streamlines when they 
cross the line source at y'—O. For small values of 7/m this 
discontinuity in vertical velocity agrees with the boundary 
condition for a simple wedge shape having a small wedge 
angle. (See fig. 2.) 

If the source strength 7 is held constant and m is allowed 
to approach zero, the wedge _ angle ultimately becomes large. 
At m=0 the line source actually satisfies the boundary con- 
dition for the circular cone (reference 7), but it is found that 
the slope of the conical boundary does not agree with the 
slope of the streamlines near 2=0 and hence the theory no 

Tfl 

longer holds. The condition y->0 thus represents the transi- 
tion from an oblique airfoil to a body of revolution and will 
be avoided in the present analysis by restricting the formulas 
to flat bodies, that is, airfoils that are thin in both longitudinal 
and transverse section. 

AIRFOIL OF WEDGE SECTION 

Over the w r edge section near the plane 2 = 0 , the formula (6) 
becomes simply 

11=7 cosh" 1 j^Tj (9) 


and is conveniently represented by the variation along a line 
parallel to the X-axis. Figure 2 shows the oblique wedge- 
shape figure corresponding to a line source with m<1.0. In 
this case the pressure field is confined to the interior of the 
Mach cone 3?—y*— £=x ,% — y' 2 — s /s =0 and the theory, 
unlike the Ackeret theory, indicates a stagnation point along 
the leading edge. (Actually, of course, the thin-airfoil 
theory shows an infinite velocity at such points, but this is 
to be interpreted as a velocity of the order of magnitude of 
the flight velocity V. The pressure to be expected along the 
leading edge is the stagnation pressure corresponding to the 
transverse velocity component.) 

Given the wedge angle measured in downstream 

sections, the source strength must vary with m according to 




m 


dz 


(ID 


(from equation (7)). Then 


A p 2 dz m x' 

r^5^ oos!l ' fi 


( 12 ) 


If m exceeds 1.0, the leading edge of the airfoil will lie out- 
side the Mach cone. In this case 


m 


Aj?_2 dz 

S ~vlx 


cos 


-i 




(13) 
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In the region between the leading edge and the Mach cone 
cos -1 i^-j is constant and equal to ; hence the pressure in 
this region is constant, that is, 


Ap 0 dz m 

y lm 2 — 1 


(14) 


Figure 3 illustrates this result. 



If m— » » a semi-infinite airfoil with its leading edge at 
right angles to the direction of flight is obtained; here 


x-my _ -y 

V(j/-ma;) ! -Kl-mV ■> ' 

and = 2 ^ wherever ^x 1 —z 2 . This value agrees with 

the Ackeret theory. 

AIRFOILS BOUNDED BY PLANE SURFACES 


The distribution of pressure over symmetrical airfoils 
bounded by plane surfaces can be obtained by superimposing 
the pressure fields for several line sources and sinks. This 
superposition is greatly simplified by the conical form of the 
pressure field for each single line source. Because of this 
form, the whole distribution in the plane z= 0 is, in effect, 
represented by a single curve. If the velocity field for a fine 
source beginning at the origin (equation (6)) is denoted by u 


and that beginning at x= — 1 is denoted by tt_i, and so forth, 
the sum 

«_!—«+! 

represents the velocity over a plate of uniform thickness 
having a beveled leading edge of constant width. (See fig. 4.) 
Similarly 

U-i—2 u+u + i 

represents the pressure field for an airfoil having diamond- 
shape cross_sections. 

The superposition required for several sources or sinks 
can be accomplished by manipulation of a single curve if it 
is remembered that u is a function of the ratio xjy. Figure 4 
illustrates this process for a source and a sink. In terms of 
the ratio x/y the separation of source and sink and honco 
the scale of the chord length continually diminishes with 
increasing distance from the root section. 

At large distances from the vertex (x'-» <» ) the expression 
(for m<1.0) 

^-^oc cosh- 1 il^-cosh- 1 j|^j (16) 
is found to approach the value 


HvSHMs) <"> 

where Q 0 Is' the Legendre function. (See reference 7.) 



Fioubi i .— Superposition of source and sink to obtain plate with beveled odge. 
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In the thin-airfoil theory for subsonic speeds it 
can be shown that if 


wcc P a (x) 




ocF^ 

dx 

(18) 

then 






uazQ n {x) 

(19) 

since Neumann’s 

I 

formula (reference 8, p. 116) 



Qn~2 t 

r & 

J-i (*-$) ? 

(20) 


may be interpreted as the integration for the velocity distri- 
bution due to an array of sources of strength 

w dZ=P n (£) d$ 

along the chord of the airfoil. The expression Q 0 — of equa- 
tion (17) thus represents the subsonic pressure distribution 
over the beveled edge. 1 

At the root section (y—0) only the forward source need 
be considered since-the airfoil surface is ahead of the Mach 
cone originating at the rear source. Here 

tt-.-tt-HOCCOSh- ,^0^)1 

occosh -1 — (21) 

and the pressure over the root section is thus constant, as 
given by the Ackeret theory, but is altered in magnitude 
by the obliquity. 

The oblique wing lying behind the Mach lines thus shows 
the Ackeret type of pressure distribution over the foremost 
section and a progressive change along the span from this 
distribution to the subsonic type of distribution. Since the 
subsonic type of distribution shows no pressure drag, there 
is a continuous falling off of the pressure drag with increasing 
distance from the root section. The pressure drag of the 
oblique wing thus arises chiefly on the foremost section, and 
it follows that the drag coefficient of the wing as a whole 
diminishes with increasing aspect ratio. It will be shown 
subsequently that the effect of cutting the wing off along a 
line y— Constant to produce a downstream tip causes a 
reduction of the pressure drag on the adjacent sections; and 
if the aspect ratio is sufficiently high, the pressure drag in 
the region of the downstream tip may actually be negative. 

If the wing lies ahead of the Mach lines (mj>1.0) the 
Ackeret type of pressure distribution occurs and a pressure 
drag arises over the whole length. In this case both u and w 
are constant over the beveled part at a distance from the origin. 

i Similarly if i\(4) di is taken as the ehordwise distribution of vorticity. 


u a P, (r) 
to cc <J« (z) 




Piocas S.— Variation of pressure distribution along span of sweptback wing, m—tan 30*. 


Slope , dzjdx 



Section parallel to flight velocity 



Distance from root section in //£? Chord lengths 

riOPBl 6.— Variation of drag coefficient with di3tanco from root section for sweptback wing. 

Wedge section; Af— M. 

The treatment thus far applies to semi-infinite cylindrical 
wings having root sections near the origin. A complete 
sweptback wing may be obtained by the addition of a sym- 
metrical or conjugate arrangement of source lines below the 
X-axis. Values of u for this conjugate arrangement may be 
denoted by u. Figures 2 and 3 show u for a single inclined 
source and figure 5 shows calculated pressure distributions 
at several sections along the span for a complete sweptback 
airfoil having beveled sections. The addition of the con- 
jugate source lines doubles the pressure at the root section, 
but this interference effect falls off rapidly along the span. 
It is noted that, as in figure 4, the most significant change 
in pressure distribution occurs along the expansion wave 
originating at the trailing edge of the root section. Figure 6 
shows the variation in pressure drag along the span for this 
airfoil obtained by integrating the ehordwise components of 
pressure at the different sections. 

The addition of a reversed source-sink distribution having 
its origin displaced to a point 0> (see fig. 7) will show the 
effect of cutting the wing off in a direction parallel to the 
direction of flight. It will be evident that the effect of such 
a tip is characterized by the subtraction of the curves u and 
is limited to the area lying within the Mach cone which 
originates at the tip. It is interesting to note that pressure 
distributions of the Ackeret type, except reversed in sign, 
are added near the tip; hence, cutting the tip off in this man- 
ner reduces the drag of adjacent sections. 


The first of this series of airfoils Is the camber shape carved to support a uniform load. 




2 ' 



Figuhe 7.— Addition of reversed source-sink distribution to produce tip. 


Figure 8 shows the pressure distributions over a rectangu- 
lar airfoil having a leading edge at right angles to the flow. 
In the triangular area ahead of the Mach cones originating 
at the tips the pressure is constant, as given by the Ackeret 
theory, whereas behind these Mach cones the pressure drops 
sharply. 

AIRFOIL OF BICONVEX SECTIONS 


Curved surfaces require a continuous distribution of sources 
and sinks alined with the generators of the surface. Each 
elementary source line causes an infinitesimal change in 
direction of the surface and hence the slope at any point may 
be obtained by adding up the effects of all sources ahead of 
that point. Thus 


dz_ v C 1 Vl— m 2 dl , 
dx~ V Jx, m 



cPz t Vl— m 2 dl 
m dx 


For airfoils of constant chord, m will be a constant and the 
integrations can be performed without difficulty. Thesimplest 
case is that of constant curvature, which leads to profiles 
formed from circular arcs. 

In order to obtain a biconvex profile, it is necessary to 
introduce finite sources of strength sufficient to form the 
desired angle of intersection of the arcs at the leading and 
trailing edges, together with a uniform distribution of sinks 
along the chord line between the two sources. These pro- 
files thus require a uniform distribution of sources or sinks, 
which may be obtained by integrating the elementary solu- 
tion for the line source (equation (6) ) , The resulting solution 

may be denoted by ^ u and is, for m< 1, 
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Floras 8.— Pressure distribution over wing or biconvex section. — — 

Q r 

This superposition may be accomplished conveniently by 

transposing and adding the typical curves u and 
fig. 9.) 

It will be found that if m is less than 1.0 the velocity 
distribution approaches, with increasing distance from the 
root section, the form given by the subsonic-flow theory for 
an airfoil of biconvex section, that is, 




At the root section, however, the form is simply that given by 
the Ackeret theory for a straight airfoil although the values 

7YI 1 

are reduced in magnitude by the factor ^=== cosh -1 — ■ 

The pressure distribution and the variation of drag along 
the span for the bilaterally symmetrical wing are shown in 
figures 10 and 11. 



Fiarag II.— Variation of drag coefficient with distance from root secHon for wings of biconvex 
section. JWf=-1.4; ^—10 percent. 


The surface obtained is one generated by a line passing 
through the fixed point Xq, y 0 and hence is a conical surface. 

The pressure over the tapered airfoil requires the integra- 
tion of 


CONICAL SURFACES 

For tapered airfoils both m and I will be functions of £. 
It is easily seen that closed surfaces can be obtained only if 
the relation between m and | is such that the line sources 
have a common point of intersection, as in figure 7. If tbia 
point is denoted by a®, y 0 



where & is the location of the vertex of the airfoil and 
dl V m d?z 
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In conclusion it should be noted that the pressures have 
been derived for an airfoil transformed according to equa- 
tions (2). The pressures at corresponding points of the 
original airfoil are to be obtained by dividing by M— 1. 

Langley Memorial Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., May £3, 1946. 


APPENDIX 

EVALUATION OF INTEGRAL OF EQUATION (8) 

For m<A.O the disturbance is zero outside the Mach cone 
and the range of integration should be extended only from 
*,= V?+P to x, that is, 




_ cosh -1 

w; v 


fr+^ dx 


(for unit source strength) . Furthermore, 


dx 


b f* . f 1 

_ J VFR dz ' 

since the integrand is zero at the lower limit. Now 


x' — x’z'sjl— m* 

bz C0Sh Vi/ ,s + 2 7 ' j ~ G/ ,! + z' 2 ) ^x n -y ri -z'' 1 


and hence the integral 


J 'x -x'z'-y/l—m 2 

,V+? dx 


(Al) 


(A2) 


(A3) 


(A4) 


must be evaluated. 

First it is noted that the integral vanishes with z except 
in the neighborhood of the Mach cone (Vz' 2 — y' 2 — z' J== 0) 
and in the neighborhood of the line source (y'=0). Near 
the Mach cone x' 2 , so that 

f —x'z'dx f —zdx e k z\ 

J (y' i +2' i )^x' 2 -y' i -z'^J 1 ] 

Since the latter integral approaches zero with z, there is no 
contribution to equation (A4) in the region of the Mach 
cone. On the other hand, near the line source y'—>0 and 
sjx' 2 —y n — z' 2 -> x' ; hence, as z'— >0, 

c -x'z'dx . f -s' , 

J (y' 2 +z' 2 )kjx' 2 -y' 3 -z' 2 J y' 2 ~z' 2 

tan -1 ^7+ Constant (A6) 


m 


The value of the integral changes from 0 to t in crossing 
over the line source at y' = 0 and is positive or negative de- 
pending on whether z' approaches zero from, the positive or 
negative side of the zy-plane. Hence 

^l—m 2 

m 


If m is greater than 1.0, 


*‘ = 00 S "vFf? s ~ (A7) 


and the flow disturbance extends outside the Mach cone to 
a region bounded by plane waves extending from the line 
source and tangent to the Mach cone. (See fig. 12.) The 
equation of these planes can be easily shown to bo y'*4-z'*=0; 
hence for m>1.0 the lower limit of integration is given by 

y'*+z'*=0 



Figure 12.— Information pertinent to evaluation of equation (8) for m>1.0. (See append!*.) 



THIN OBLIQUE AIRFOILS AT SUPERSONIC SPEED 


275 


or 

Xi — 

m 

(AS) 

Then 




(A9) 


In this case u does not go to zero at the lower limit but is 
equal to r. In all other regions, however, the integral 
approaches zero uniformly with 2 as in the preceding case; 
hence 



u dx 



= v^l (A10) 

as before. 
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